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Supplementary Text 

Explorations of additional candidate genes associated with root length 

A novel association between mean root length and At5g45120 

We identified an SNV on chromosome 5 (18,241,640 bp - SNP 5_18241640) in the epistatic 
GWAS analysis through its interaction with SNV 3_66596 (Table 2). Further analyses showed 
that this SNV was in LD with two missense SNVs and four synonymous SNVs in the gene 
At5g45120. The insertion mutant of At5g45120, however, showed wild-type root length (Figure 
4), possibly due to the limitations of single-mutant analysis in validating epistatic interactions. 

At5g45120 encodes a member of the little studied eukaryotic aspartyl protease family (Table 4) 
[1]. The two missense variants, T>T and D->A, were positioned 5 bp apart and were either 
adjacent to or within a 'phosphorylation hotspot', suggesting that post-translational 
modifications may be altered in some accessions [2,3]. Although mutating At5g45120 in a Col-0 
genetic background did not affect mean root length, future exploration of this mutant or the 
minor allele in alternative genetic backgrounds may reveal whether At5g45120 affects root 
development. 

At3g01180 and At3g01 185 are positional candidate genes for mean root length 
The SNV 3_66596 was associated with mean root length in the epistatic GWAS through 
interactions with i) two SNVs at the ATL5 locus (Table 2) and ii) a SNV on chromosome 5 
(5_18241640). The SNV resides in the gene At3g01185, a 148 amino acid protein of unknown 
function (Table 4) [4], where it results in a synonymous change. This SNV also resides in the 
adjacent promoter region of the gene At3g01 180, which is highly accessible by DNase I [5]. As 



this SNV is not in LD with any other SNVs in the 1001 Genomes re-sequencing data, it might be 
the causal SNV for the observed association, either by directly affecting the function of 
At3g01 185 or by altering the expression of At3g01 180. 

We phenotyped the roots of an insertion mutant for At3g01 185, but failed to find significant 
effects on root growth (Figure 4), implying that disruption of this gene alone is not sufficient to 
affect root development. Unfortunately, an insertion mutant for At3g01 180 was not available. 
Further studies, including generating and phenotyping At3g01 180 mutants, will be required to 
either validate or exclude one of these candidate genes. 

Epistatic interaction with a transposon-rich region on chromosome 1 

SNV 1_17257526 interacted with the NAC6 locus on chromosome 5. This SNV is located in the 
gypsy-like transposon Atlg46624 and was not in LD with any other SNVs. The 20kb region 
surrounding the peak association encodes 16 transposable elements and many regulatory 
elements as identified by DNasel-accessible chromatin [5]. How a polymorphism in this 
complex and potentially volatile region would affect root length by itself, or through epistatic 
interactions, is unclear and deserves further investigation. 
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Figure SI. Experimental design to increase accuracy of variance estimates. For 93 

accessions, a total of 210 seedlings were grown in three replicates of 70 individuals. For each 
replicate, the mean and variance was estimated. The mean and variance ( 0*0= ) of the variance 
across three replicates was also calculated. 



accession 


stock 
number 


Ag-0 


CS22630 


An-1 


CS22626 


Bay-0 


CS22633 


Bil-5 


CS22578 


Bor-1 


CS22590 


Bor-4 


CS22591 


Br-0 


CS22628 


Bur-0 


CS22656 


C24 


CS22620 


CIBC-17 


CS22603 


CIBC-5 


CS22602 


Col-0 


CS22625 


N13 


CS22491 


Ct-1 


CS22639 


CVI-0 


CS22614 


Eden-1 


CS22572 


Eden-2 


CS22573 


Edi-0 


CS22657 


Ei-2 


CS22616 


Est-1 


CS22629 


Fab-2 


CS22576 


Fab-4 


CS22577 


Fei-0 


CS22645 


Ga-0 


CS22634 


Got-22 


CS22609 


Got-7 


CS22608 


Gu-0 


CS22617 


Gy-0 


CS22631 


HR-10 


CS22597 


HR-5 


CS22596 


Kas-2 


CS6751 



accession 


stock 
number 


Kin-0 


CS22654 


Knox-10 


CS2241 1 


Knox-18 


CS22567 


Kondara 


CS22651 


Kz-1 


CS22606 


Kz-9 


CS22607 


Ler-1 


CS22618 


LL-0 


CS22650 


L6v-1 


CS22574 


L6v-5 


CS22575 


Lp2-2 


CS22594 


Lp2-6 


CS22595 


Lz-0 


CS2261 5 


Mr-0 


CS22640 


Mrk-0 


CS22635 


Ms-0 


CS22655 


Mt-0 


CS22642 


Mz-0 


CS22636 


Nd-1 


CS2261 9 


NFA-10 


CS22599 


NFA-8 


CS22598 


Nok-3 


CS22643 


OM62-1 


CS22584 


OM62-3 


CS22585 


Oy-0 


CS22658 


Pna-17 


CS22570 


Pro-0 


CS22649 


Pu2-23 


CS22593 


Pu2-7 


CS22592 


Ra-0 


CS22632 


Ren-1 


CS2261 0 



accession 


stock 
number 


Ren-1 1 


CS2261 1 


Rmx-A02 


CS22568 


Rmx-A180 


CS22569 


RRS-10 


CS22565 


RRS-7 


CS22564 


Se-0 


CS22646 


Shakdara 


CS22652 


Sorbo 


CS22653 


Spr1-2 


CS22582 


Sq-1 


CS22600 


Sq-8 


CS22601 


Tamm-2 


CS22604 


Tamm-27 


CS22605 


Ts-1 


CS22647 


Ts-5 


CS22648 


Tsu-1 


CS22641 


UII2-3 


CS22587 


UII2-5 


CS22586 


Uod-1 


CS22612 


Uod-7 


CS22613 


Van-0 


CS22627 


Var2-1 


CS22580 


Var2-6 


CS22581 


Wa-1 


CS22644 


Wei-0 


CS22622 


Ws-0 


CS22623 


Ws-2 


CS22659 


Wt-5 


CS22637 


Yo-0 


CS22624 


Zdr-1 


CS22588 


Zdr-6 


CS22589 



Table SI. A. thaliana accessions 



phenotyped for GWA. 
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Figure S2. Root length variance varies independently of root mean and can be calculated 
globally or per replicate. A) There is a strong correlation between the total variance estimated 
from -210 individuals and mean variance from three replicates of -70 individuals (p < 2.2e-16, 
R = 0.98). B) There is no significant correlation between root length mean and root length 
variance (p = 0.11, R 2 = 0.16, Pearson's). 
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Figure S3. No associations are found with root length using an additive model. A) A GWAS 
was completed under an additive model. No SNPs reached the Bonferroni threshold. B) Using a 
whole-genome generalized ridge regression model in which SNPs were random effects, no SNPs 
reached the Bonferroni threshold. 



7 -| 




Chromosome 

Figure S4. No loci are detected for root length variance using GWAS with an additive 
model. A GWAS was completed under an additive model. No SNPs reached the Bonferroni 
threshold. 
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Figure S5. The same pattern of effects is observed for all significantly interacting loci. A-D) 
Root mean is displayed for the four pairs of significantly interacting loci. The major allele is 
indicated by -1 and the minor allele is indicated by 1. The combination of two minor alleles (1,1) 
always has the largest root length compared to the other three allele combinations, which are 
more similar to one another. A) Root mean for allele combinations atl_17257526 and 
5_15862026 /5_15862525. B) Root mean for allele combinations at 3_66596 and 
3_9272294/3_9273674. C) Root mean for allele combination at 3_66596 and 5_18241640. D) 
Root mean for allele combination at 3_10891 195 and 5_1027939. 
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T-DNA line 


Insertion 
location 


At3g01185 


SALK 


_033165 


promoter 


At3g01200 


SALK 


_035653 


Exon 1 of 3 


At3g25520 


SALK 


_1 35037 


5'UTR 


At3g25530 


SALK 


_05741 0 


Intron 3 of 7 


At3g25540 


SALK 


_073240 


Exon 1 of 5 


At5g45120 


SALK 


_030666 


5'UTR 


At5g39610 


SALK 


_090154 


Exon 3 of 3 


At5g39620 


SALK 


_096950 


Intergenic 


At3g28865 


SALK 


_089352 


Exon 1 of 1 


At3g28865 


SALK 


_009759 


Exon 1 of 1 


At3g28880 


SALK 


_070443 


Exon 17 of 19 


At5g 03840 


SALK 


_1 42051 


5' UTR 


At5g03850 


SALK 


_1 12775 


Exon 1 of 1 



Table S2. T-DNA lines tested for root phenotypes. 



